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Abstract 


Homogeneous and isotropic, nonsingular, bouncing world models are designed to evade the initial 
singularity at the beginning of the cosmic expansion. Here, we study the thermodynamics of the 
subset of these models governed by general relativity. Considering the entropy of matter and 
radiation and considering the entropy of the apparent horizon to be proportional to its area, we 
argue that these models do not respect the generalized second law of thermodynamics, also away 
from the bounce. 
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I. INTRODUCTION 


The well-know, and partially successful, cosmological standard model, based on Einstein 
gravity, assumes a homogeneous and isotropic universe whose energy density is supplied by 
matter and fields that satisfy the strong energy condition |l|. The popular Einstein and 
de Sitter cosmology is a case in point. In this class of models an inevitable singularity 
occurs at the beginning of the expansion. As f —)■ 0 the energy density grows beyond 
whatever bound, world lines cannot be continued any longer and, worse of all, the spacetime 
description breaks down. This holds also true if at very early times, before the conventional 
radiation dominated phase, the expansion of the universe is driven by some held whose 
energy density, essentially constant, overwhelms every other form of energy j^. 


To evade this uncomfortable, philosophical, implication three main possibilities were ad¬ 
vanced: (i) The initial singularity arises chiehy from ignoring quantum effects which are to 
prevail at some point before the singularity can occur . (ii) The universe began expanding 
from a hnite initial size safely larger than Planck’s length (thus avoiding as well the quantum 
dominated regime). In this scenario, known as the “emergent” universe, an Einstein-type 
static phase replaces the initial singularity —see e.g. js, 6|. (iii) Nonsingular bouncing 
universes. According to this proposal the universe experiences either just one contracting 
phase followed by an expanding one, or a sequence of expanding and contracting phases 
connected, in either case, to each other by a bounce such that the scale factor, aft), of the 
Friedmann-Robertson-Walker (FRW) metric neither vanishes nor diverges a -Q — see 


and 


14| for reviews. Since the singularity theorems of Penrose and Hawking Q forbid 


bounces to originate from normal matter the physics behind them is not straightforward 
and requires the presence of some source of energy that violates the null energy condition. 
It is therefore understandable that the bounce mechanism remains a subject of debate as 
yet 91] . Nevertheless, nonsingular bouncing may arise in de Sitter-like models with closed 
spatial sections as well as in supergravity-based models [^. Figure H] shows the evolu¬ 
tion of the scale factor and Hubble factor, H = a/a —in arbitrary units—of a nonsingular, 
bouncing toy model. In spite of its simplicity, it encapsulates the salient kinematics features 
of this class of models. 

Given the strong connection between gravity and thermodynamics 18j-2l||, for a cosmo- 
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Figure 1. Qualitative evolution of the scale factor, solid (black) line, and the Hubble function, 
dashed (red) line, in arbitrary units, of a nonsingular bouncing FRW universe. Here, o = 2 + sin(t). 


logical model that fits reasonably well the observational data it will be a bonus to comply 
with the laws of thermodynamics, especially the second law. The latter formalizes our 
daily experience that macroscopic systems, left to themselves, spontaneously tend to ther¬ 
modynamic equilibrium or remain in it, if undisturbed. According to it, the entropy S of 
isolated systems never decreases. S' > 0, and must be concave at least in the last stage 
of approaching equilibrium 2^. In the mid-1970s this was proved for gravity dominated 


systems, namely, black holes plus their environment [23 


24| . and later extended to cosmic 


scenarios consisting in a de Sitter event horizon plus the fluid enclosed by it 2^, l26 |. 


In the case of homogeneous and isotropic universes, the entropy measured by a comoving 
observer splits into two parts. S' = -b S'/. The hrst one is the entropy of the apparent 
horizon 271 

S^ = kBA/{Ai%) , (1) 


where A = denotes its area and 


fA = 1/VH^ + ka-^ 


( 2 ) 


denotes the horizon radius 


27l | . This is in keeping with the fact that the entropy of systems 


dominated by gravity does not vary as its volume but as the area of the surface that bounds 
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this volumq^. The second part, Sf, corresponds to the entropy of fluids (e.g. matter and 
radiation) enclosed by the horizon. In many instances of interest the latter is negligible as 
compared with the former; this is only natural given the smallness of the factor. (We 
do not consider the entropy of scalar helds as we assume these fields are in a pure quantum 
state and consequently their entropy vanishes altogether). In the particular case of the 
present universe the entropy of the horizon is larger than that provided by supermassive 


black holes, stellar black ho 
of magnitude, respectively 


es, relic neutrinos and CMB photons by 18, 25, 33 and 33 orders 


At this point, we feel it expedient to summarize the derivation of the radius of the apparent 
horizon, Eq. (|5]). We begin by rewriting the FRW metric, 

ds^ = —dt^ + a?{t) —+ a?{t) dVt^ (3) 

1 — kr‘^ 

(where k = +1, 0, —1, denotes the spatial curvature index and hi is the unit two-sphere), as 


ds^ = hahdx^ dx’’ + r‘^{x) dhl^ . 


( 4 ) 


Here, x^ = t, x^ = r, and hab = diag 


— 1 

’ l—kr'^ 


The radius of the dynamical apparent 


horizon is set by the condition h°‘^dafdbf = 0. A straightforward calculation produces Eq. 
(!2|) . To better understand the meaning of the horizon note that the expansion of the ingoing 
and outgoing null geodesic congruences are given by 


1 I 1 I 

9in = H - - \ 1 -- and 9out = H + - \1 -- , (5) 

r \ r \ 

respectively. A spherically symmetric spacetime region will be called a “trapped” region if 
the expansion of ingoing and outgoing null geodesics, normal to the spatial two-sphere of 
radius r centered at the origin, is negative. By contrast, the region will be called “anti¬ 
trapped” if the expansion of both kind of geodesics is positive. In normal regions outgoing 
null rays have positive expansion and ingoing null rays, negative expansion. Thus, the anti¬ 
trapped region is given by the condition f > + ka~‘^)~^^‘^. Clearly, the surface of the 

apparent horizon is nothing but the boundary hypersurface of the spacetime anti-trapped 
region. Therefore, it is only natural to associate an entropy to the apparent horizon since 


^ This is, e.g., the case of the entropy of black holes. 
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one can regard the latter as a measure of one’s lack of knowledge about the rest of the 
universe. Further details concerning the entropy of the apparent horizon can be found in 
j^ . Obviously, in the case of an exact de Sitter expansion, the apparent and event horizons 
coincide. 


Models falling into the emergent class fulhll the generalized second law (GSL) 29|, i.e., 
Sa + Sf > 0, as well as the concordance ACDM model |30|, IM]- Likewise, the models of 
references |32| and [33[ satisfy the said law [3^. Both scenarios evade the big bang singu¬ 


larity and evolve from an initia’ 
the far future. In the hrst one. 


de Sitter expansion to a hnal, also de Sitter, expansion in 


32| , a gravitationally induced production of particles drives 


the current cosmic acceleration. In the second one. 


33| . the latter arises from a dynamical 


cosmological “constant” that depends on the Hubble factor in a manner dictated by the 


renormalization group j^| 
Einstein and de Sitter [2| 
fail to satisfy this law [30|, 


By contrast, ever expanding matter-dominated universes (as the 
and phantom-dominated cosmologies obeying Einstein gravity 


The target of this paper is to explore whether the GSL is satished by the nonsingular 
bouncing class of world models that obey general relativity, in their evolution outside the 
bounce itself. As mentioned above, during the bounce the null energy condition, p + p > 0, 
is violated; consequently the second law is expected to be violated as well. It is therefore 
worthwhile to ask whether the transgression of the second law is conhned just to the bounce 
or it also occurs in other phases of the universe evolution. This is the question we aim to 
answer. 


Next section focuses on a homogeneous and isotropic cosmological model, sourcered by a 
massive, complex, scalar held, that display a nonsingunlar bounce [^. Section III studies 
from an overall viewpoint the class of models that incorporate matter and radiation as 
sources of the gravitational held though, regrettably, these models do not have analytical 
solutions. Section IV presents our conclusions and hnal comments. 
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II. NONSINGULAR STAROBINSKY’S MODEL 


Let us consider a model proposed by Starobinsky whose only energy component is a ho 


mogeneous and massive, complex, scalar held [^. This is the simplest model featuring a 
nonsingular bounce that has analytical solutions for the scale factor with a nonzero measure 
in the space of initial conditions —in this last respect; see also 


The model obeys Friedmann’s and energy conservation equations 


k SttG 


^ = 


{\^\^ + , 


ip + 3Hip + m'^ip = 0, 

where m is the quantum mass of the held and fc = +1 in order for Eq. 


( 6 ) 

( 7 ) 

to allow a bounce. 


At the epoch of maximum expansion (—Omax -C t —m~^) the universe contracts according 
to a{t) = where A is a positive-dehnite constant. In this era the variation with time 

of the area of the apparent horizon is 


A 


—dvr 



+ A-2|t|-4/3)2 


< 0 . 


( 8 ) 


For t S> m ^ but t < m ^(2 InmOmax)^'^^ the scale factor obeys a{t) 
C = im?lQ j^. Consequently, 


A 


-dvr t 


on2 I 4C'exp(2Ct2) 


< 0 . 


oi exp(—with 


( 9 ) 


Figure [2] shows the evolution of the scale factor (though, in reality, the solution is valid only 
for a portion of the graph) as well as the evolution of the entropy of the apparent horizon. 
As we see, the entropy of the horizon diminishes, at least, in some interval of the contracting 
stage. Since neither matter nor radiation sources the gravitational held the total entropy 
decreases in such periods. 


III. MODELS WITH MATTER AND RADIATION 

We next consider nonsingular bouncing models that include matter and radiation. Since we 
are not aware of analytical solutions for the scale factor our study will necessarily be from 
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Figure 2. Evolution of the scale factor —solid (black) line —given by a{t) = ai exp(—and 
the entropy of the horizon —dashed (blue) line —for t > 0. 


a general viewpoint. In the following we adopt the toy model of Fig. [3l 


Inspection of the said hgure reveals that irrespective of the sign of the curvature index, k, 
of the FRW metric (ED, the entropy of the apparent horizon S'_ 4 , oscillates so that it presents 
alternated phases of increase and decrease. To fulhll the GSL, the entropy of matter and 
helds inside the horizon should compensate for the stages of decrease. Even though this 
can hardly be the case this may not be necessarily true near the turning points (where a 
vanishes). In these points, it may depend on the details of the specihc bouncing model. So, 
except possibly for the turning points, the entropy of matter and helds inside the horizon, 
most likely, will not compensate. 


In Fig. |3]it is seen that when the scale factor goes from a maximum to the next minimum, 
the area, A = + ka~‘^) (or what amounts to the same, the entropy), of the k = +1 

apparent horizon decreases. This is also true for the fc = 0 apparent horizon when the scale 
factor goes from a maximum to the consecutive inhection point. Next we shall check if 
matter and radiation may compensate for the decrease in entropy. 


Note that, as mentioned above, in most instances, the entropy of the horizon overwhelms any other form 


of entropy 


2a. 
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Figure 3. Evolution of the scale factor —solid (black) line —and the horizon entropy, in arbitrary 
units, for A: = 0 —dashed (blue) lines —, and A: = +1 —dotted (red) line. Here, a = 2 + sin(A). 

The entropy of pressureless matter, as well as of massless radiation, inside the ap parent 
horizon varies as the number of particles in there, i.e., Sm = |37|. Since 

Ni = (47r/3) rij r3i, where the number densities obey Um = nmoci~^ and = n^Qa~‘^ for 
matter and radiation, respective!}^, one obtains 


Sm ^mO ( 

a 


a a. 


and 


= —ATiksn^o ( — 

CL 


H 


rAV 1 

— «+ TT 

a J 6a 


( 10 ) 


( 11 ) 


Inspection of the right-hand-side of the last two equations and Fig. [3l reveals that the en¬ 
tropy of the fluids (matter and radiation) enclosed by the apparent horizon also diminishes 
in the scale factor interval between each maximum and the subsequent inflection point (note 
that there both d and d are negative quantities). Thereby, in the said interval the total 
entropy (that of the horizon plus that of relativistic and non-relativistic particles inside the 
horizon) decreases, Sa + Sm + S'.y < 0. 

In addition, as is well known, when matter evolves from an uniform distribution to a 
clumped one to form bound objects (stars, galaxies, and so on) its entropy also diminishes 


^ For simplicity we have set oq = 1. 





















since spatial matter distribution goes from a disordered state to a less disordered one. In 
expanding phases of an oscillating universe this decrease might be counterbalanced by the 
entropy gain associated to the volume increase of the universe, but not in contracting phases. 


Likewise, if one includes the “entropy” associated to the perturbed FRW universe, recently 


introduced by Clifton et ah 


, based on the Bel-Robinson tensor, the situation does 


not change since the said entropy is proportional to (see Eq. (54) in 38|]). Then in 
contracting phases this entropy diminishes as well. 


For a bounce (or a sequence of bounces) to occur in the class of models considered in this 
section, an energy component that violates the null energy condition (NEC) must also be 
present. Thus, one may wonder if this component may restore the GSL in this context. 
Typically, helds that break the NEC are scalar helds of phantom type. As mentioned in the 
Introduction, it is generally assumed that scalar helds have vanishing entropy. Upon this 
usually accepted assumption, the presence of this held will not modify our results in any 
way. However, one may still argue that the said held might, in reality, be an ehective one 
that corresponds to a mixture of helds, every single one being in a pure quantum state. In 
such a case, this ehective phantom held should be entitled to an entropy. Now, the entropy 
of any system is an increasing function of its energy and volume. In the contracting phase, 
at variance with matter and helds that fulhll the NEC, the energy density of the phantom 
held decreases. Since the volume of the apparent horizon diminishes as well, the entropy of 
this ehective, phantom held, decreases too. Once again, our result stands. 


Altogether, our analysis suggests that the GSL is not satished by homogeneous and isotropic, 
nonsingular, bouncing models obeying Einstein relativity. 


IV. CONCLUDING REMARKS 

Assuming that the entropy of the universe is proportional to the area enclosing the dy¬ 
namical apparent horizon, and taking into account the contributions by the matter and 
radiation, we investigated if nonsingular bouncing models satisfy the generalized second 
law of thermodynamics. We found that homogeneous and isotropic nonsingular bouncing 
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models, governed by general relativity, do not comply with the GSL also away from the 
bounce. This is not only an interesting result in itself, it also begs the question as to whether 
more elaborated bouncing models respect this law. 

One may wonder whether our results would change if we have used the event horizon rather 
than the apparent horizon. The answer is: no. The entropy of the former horizon is also 
given by his area, Aeh where the radius of the horizon reads leh{t) = a{t) /“dt/a(t). 

Consequently, in contracting phases this entropy decreases as well. 


Similarly, the “entropy of the gravitational field” introduced four decades ago by Penrose 


39, 


40l |. based on Weyl’s tensor does not change our results because, in the case of non- 


perturbed FRW universes, the Weyl tensor vanishes identically and so does this entropy. At 
any rate, it faces well known difficulties; see e.g. 


41 


42|. 


One might still argue that the failure to fulhll the GSL arises from the fact that we are 


observing just one universe and that we did not 
from the principle of indifference of Laplace 


ake into account cosmic variance. However, 


431, it is hard to see how cosmic variance will 
alter this result for the subset of the multiversqj obeying general relativity. In any case, this 
depends on how one would choose a measure in the space of solutions, a subject that lies 
beyond the scope of this work. 


Finally, bouncing models based on gravity theories that go beyond general relativity, or that 


introduce extra dimensions (see [L 


posals by Steinhardt et ah [44 


IQ 


and references therein, as well as the recent pro- 


45| and Ijjas and Steinhardt 


46|j. should also be considered 


from the thermodynamic standpoint. However, for the time being, it is rather unclear how 
to define a meaningful entropy in such models. 


We conclude by saying that, while the aforesaid connection between gravity and thermo¬ 


dynamics. 


18| - 21| . does not fully ensure us that the universe behaves as a thermodynamic 


system (i.e., one that obeys the thermodynamic laws), given two models that £t equally well 
the observational data, one respecting the GSL while the other does not, the former will 


By multiverse we mean the set of all model universes. 
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likely be preferred. 
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